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Popular summary. 

Aerosols, tiny solid or liquid particles suspended in the atmosphere, were once only a 
side note in the Atmospheric Sciences. Today we realize the importance of aerosols in 
instigating or mitigating climate change, in modifying clouds and large-scale 
precipitation patterns and in affecting human health. Unlike greenhouse gases, which are 
well-mixed and long-lasting in the atmosphere, aerosols are temporally and spatially 
variable with lifetimes of a few days to a few weeks. Their transient natures make 
aerosols difficult to characterize and their effects on climate, hydrology and health 
difficult to model. Satellites provide the best means to observe the global aerosol system 
and narrow the uncertainties associated with aerosol characterization, but the satellite 
observations must be sufficiently accurate to be useful. The MODerate resolution 
Imaging Spectroradiometer (MODIS) aboard both NASA's Terra and Aqua satellites 
provides a unique tool to discern the global impact of aerosols. The products derived 
from MODIS data include aerosol optical thickness, which is a measure of aerosol 
amount, as well as products that describe the size of the aerosol particles. The MODIS 
aerosol retrievals are continuously evaluated against ground-truth of an existing global 
network of highly accurate instruments (AERONET). The results show an accuracy for 
the MODIS aerosol products that will sufficiently narrow the uncertainty of global 
aerosol characterization. Furthermore, the MODIS derivation of aerosol particle size aids 
in discriminating between man-made aerosol and naturally produced aerosols. This is a 
major step forward in narrowing the uncertainties associated with estimating the total 
anthropogenic effect on climate. 
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Abstract. The MODerate resolution Imaging Spectroradiometer (MODIS) aboard both NASA’s 
Terra and Aqua satellites is making near global daily observations of the earth in a wide spectral 
range. These measurements are used to derive spectral aerosol optical thickness and aerosol size 
parameters over both land and ocean. The aerosol products available over land include aerosol 
optical thickness at three visible wavelengths, a measure of the fraction of aerosol optical 
thickness attributed to the fine mode and several derived parameters including reflected spectral 
solar flux at top of atmosphere. Over ocean, the aerosol optical thickness is provided in seven 
wavelengths from 0.47 pm to 2.13 pm. In addition, quantitative aerosol size information includes 
effective radius of the aerosol and quantitative fraction of optical thickness attributed to the fine 
mode. Spectral aerosol flux, mass concentration and number of cloud condensation nuclei round 
out the list of available aerosol products over the ocean. The spectral optical thickness and 
effective radius of the aerosol over the ocean are validated by comparison with two years of 
AERONET data gleaned from 133 AERONET stations. 8000 MODIS aerosol retrievals co- 
located with AERONET measurements confirm that one-standard deviation of MODIS optical 
thickness retrievals fall within the predicted uncertainty of At=±0.03±0.05t over ocean and 
At=±0.05±0.15t over land. 271 MODIS aerosol retrievals co-located with AERONET 
inversions at island and coastal sites suggest that one-standard deviation of MODIS effective 
radius retrievals falls within Ar eff = ±0.11 pm. The accuracy of the MODIS retrievals suggests 
that the product can be used to help narrow the uncertainties associated with aerosol radiative 


forcing of global climate. 
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1. Introduction 

The Chesapeake Lighthouse and Aircraft Measurements for Satellites (CLAMS) field 
experiment was designed to aid the development and evaluation of satellite algorithms that retrieve 
geophysical parameters important to the earth's radiative balance and estimates of global change. 
Aerosols are one of those important geophysical parameters that are active players in the earth's energy 
balance and hydrological cycle. These suspended airborne particles scatter solar radiation back to 
space, absorb solar radiation in the atmosphere and shade the earth's surface. Airborne particles act as 
cloud condensation nuclei, entering into cloud processes and thereby change cloud reflectivity and the 
hydrological cycle (Twomey, 1977; Rosenfeld and Lensky, 1998). Aerosols also affect human health 
and reduce visibility (Samet et al. 2000). Some aerosol types are natural such as wind blown desert 
dust or sea salt caused by breaking waves. Other aerosol types are created from human activities such 

as urban/industrial pollution and biomass burning. Unlike C0 2 , another atmospheric pollutant input 

into the atmosphere from human activity, aerosols are not well mixed in the atmosphere, and because 
of their spatial and temporal variability, the uncertainty of estimating human-induced aerosol forcing 
on climate and the hydrological cycle is huge. Therefore, characterizing global aerosol distribution 
presents one of our major challenges today (Kaufman, et al. 2002). 

Operational production of aerosol information from long-term satellites provides a means to 
achieve a global and seasonal characterization of aerosol. Satellites view the entire earth and produce a 
global image, thus resolving the spatial patterns resulting from the spatial inhomogeneities of aerosol 
sources. Daily global images from polar orbiting satellites (Husar et al., 1997; Torres et al. 2002) and 
more frequent imagery from geostationary satellites (Prins et al. 1998) resolve the temporal patterns 
resulting from the short lifetimes of aerosols, which are on the order of a few days to a week. 
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The MODerate resolution Imaging Spectroradiometer (MODIS) is a new sensor with the ability 
to characterize the spatial and temporal characteristics of the global aerosol field. Launched aboard 
NASA's Terra and Aqua satellites in December 1999 and May 2002, MODIS has 36 channels spanning 
the spectral range from 0.44 to 15 pm. The aerosol retrieval makes use of seven of these channels 
(0.47 - 2.13 pm) to retrieve aerosol characteristics, and uses additional wavelengths in other parts of 
the spectrum to identify clouds and river sediments (Ackerman et al., 1998 ; Gao et al., 2002; Li et al. 
2003) . Cloud identification is also aided by the 500 m spatial resolution (Martins et al. 2002). Because 
of the wide spectral range unavailable on previous satellites, MODIS has the unique ability to not only 
retrieve aerosol optical thickness with greater accuracy, but also to retrieve parameters characterizing 
aerosol size (Tanre et al. 1996; Tanre et ah 1997). The results section of this paper show that MODIS's 
ability to separate aerosols by size can be used as a proxy for separating human-generated aerosol from 
natural sources, which aids substantially in estimating global human-induced aerosol forcing 
(Kaufman, et ah 2002). 

The first MODIS instrument was launched aboard the Terra satellite at the end of 1999 and 
began transmitting data at the end of February 2000. Algorithms were in place, designed to use the 
observed radiances to derive many important aerosol products. Early comparisons of the retrieved 
aerosol parameters with ground-based validation data, showed remarkable agreement between the two 
types of data (Chu et ah 2002; Remer et ah, 2002), but also showed us situations in which the 
algorithms could be improved. Almost immediately, the algorithms were modified to reflect a better 
understanding of the instrument's capabilities and the nature of aerosols and clouds. In a companion 
study in this special issue the MODIS aerosol algorithm over ocean is compared with an independent 
aerosol retrieval algorithm applied to the same data set of MODIS radiances (Ignatov et al., this issue). 
In another companion paper in this special issue, the MODIS retrievals over land and ocean are 
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evaluated regionally for the specific t im e and location of the CLAMS field study using the additional 
resources available during the CLAMS intensive observing period (Levy et al., this issue). However, 
in the present study we take a global view. We give a comprehensive description of the MODIS 
aerosol algorithms, highlighting the changes that were implemented post-launch. We describe the 
wealth of aerosol products derived from MODIS data and available to any user. Lastly, we show some 
of the global comparisons to ground-based data as validation for the products previously described. 

2. MODIS Aerosol Algorithms 

The MODIS aerosol algorithm is actually two entirely independent algorithms, one for deriving 
aerosols over land and the second for aerosols over ocean. Both algorithms were conceived and 
developed before Terra launch and described in depth in Kaufman, et al. (1997b) and Tanre, et al. 
(1997). In addition, Levy, et al. (2003) provide a more recent description of the over ocean retrieval 
algorithm. Both the land and ocean aerosol algorithms rely on calibrated, geolocated reflectances 
provided by the MODIS Characterization Support Team (MCST), identified as products MOD02 and 
MOD03 for Terra MODIS products and MYD02 and MYD03 for the Aqua MODIS products (MCST 
2000; MCST 2002). Ignatov et al. (this issue) provides a good discussion of these reflectances and 
possible errors associated with them. The Level lb reflectances along with the MODIS cloud mask 
product identified as MOD/MYD35 (Ackerman et al. 1998), a water vapor product identified as 
MOD/MYD07 and auxiliary data from NCEP provide the input for the algorithms. The MOD/MYD35 
cloud mask product also supplies the earth's surface information that identifies whether a pixel is a 
"land" pixel or a "water" pixel. Although the algorithm inputs MOD/MYD07 and the NCEP data, it 
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can run successfully without these supplements by using climatology for first guess water vapor and 
ozone profiles. 

The theoretical basis of the algorithms has not changed from inception, although some of the 
mechanics and details of the algorithms have evolved. MODIS data is organized by collections. A 
collection consists of products that were generated by similar, but not necessarily the same, versions of 
the algorithm. A complete history of changes to the algorithm over the course of the MODIS mission 
can be found at http:/7modis-atmos.gsfc.nasa.gov/MODQ4 L2/historv.htm l In this section we will 
leave the explanation of the theoretical basis of the algorithms to the earlier references, and instead 
focus on the mechanics of the V4.1.3 algorithm presently in operation, highlighting the changes made 
since 1 997. 

2.1 The Land Algorithm 

Figure 1 illustrates the mechanics of the land algorithm. The MODIS Level lb image, called a 
granule, consists of a 5-minute swath of MODIS -measured radiances or reflectances. These 
reflectances are corrected for water vapor, ozone and carbon dioxide before the algorithm proceeds. 
The first step in deriving aerosol products over land is to organize the measured reflectances of the 
three MODIS channels used in the procedure: p 047 , p 0 66 and p 213 . All three channels are organized into 
no min al 10 km boxes corresponding to 20 by 20 or 400 pixels for each box. This organization requires 
the 250 m resolution 0.66 pm channel to be degraded to 500 m in order to match the resolution of the 
other two channels. The 400 pixels in the box are evaluated 

pixel by pixel to determine if the standard MODIS cloud mask (MODIS product MOD35) has 
identified the pixel as cloudy, as snow/ice or as water. Since launch additional masking has been put 
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in place. Because the algorithm is sensitive to small subpixel patches of snow/ice, now all 8 pixels 
contiguous to a pixel identified as 'snow/ice' by MOD35 will also be labeled as 'snow/ice'. The pixels 
are further screened for subpixel water by determining the value of the Normalized Difference 
Vegetation Index (NDVI) for each pixel. Values of NDVI less that 0.10 are identified as containing 
subpixel water and are excluded along with cloudy and snowy pixels from the remainder of the 
algorithm. 

The second step is to select dark targets and determine their surface reflectance. This process is 
described in (Kaufman, et al. 1997b), with the following modifications. The algorithm now includes 
brighter surfaces, which expands the geographical extent of the land retrieval. The reflectance at 3.8 
pm is no longer considered, and dark pixels are selected based only on their reflectance at 2.13 pm. 
To be selected, a pixel must fall within the range of 0.01 < p 213 < 0.25. Experimentation with the 
operational retrieval showed us that p 2 13 values as bright as 0.25 gave us the same accuracy as the more 
conservative value of 0.15 initially proposed in (Kaufman, et al. 1997b)). The pixels remaining after 
masking and dark target selection are then sorted in terms of their visible reflectance, p° 66 . The pixels 
with the darkest 20% and brightest 50% of p 0 66 are discarded. The reason is to eliminate remaining 
pixels possibly contaminated by cloud shadows or odd surfaces at the dark end or residual cloud 
contamination and odd surfaces at the bright end. The possibility of residual bright cloud 
contamination is more common than cloud shadows, thus the filter is skewed towards permitting more 
dark pixels than bright ones. The remaining 30% of the pixels will be the ones used in the regular 
retrieval path, labeled Path A in Figure 1, but only if there are at least 12 of these pixels remaining 
from the original 400 in the 20 by 20 box. The mean measured reflectance is calculated from these 12 
or more dark target pixels in the three wavelengths (p 047 , p 0 ' 66 and p 2 ’ 3 ). The surface reflectances at 
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0.47 pm and 0.66 pm (p s 0 ' 47 , p s 0 ' 66 ) are derived from the mean measured p 2 13 value using the empirical 
relationships 

p s 047 =0.25 p 2 l3 ; p s 066 =0.50 p 213 (1) 

as described in Kaufman, et al., (1997b) and Kaufman, et al. (1997c). A retrieval following Path A is 
given a Quality Control value of 3, ‘very good'. An alternative path used for brighter surfaces will be 
described below. 

The third step in the standard procedure is to choose an aerosol model and make a retrieval. 
The estimated surface reflectances (p s 047 , p 0 66 ) and the measured mean top-of-atmo sphere 
reflectances (p 047 , p 0 66 ) are used as input into the Continental model Look-Up Table (LUT) to 
retrieve values for the aerosol optical thickness at 0.47 and 0.66 pm (x 047 , x 066 ). In the land algorithm, 
the two wavelengths are derived independently. The retrieved optical thicknesses along with the 
Continental model's single scattering albedoes (m 0 ° 47 , © 0 °' 66 ) and phase functions at the appropriate 
scattering angle (P 047 , P 0 66 ) are used to calculate the path radiance in each wavelength using the 
following relationship: 

p o °- 47 = co o a47 x 0 - 47 P 0 - 47 ; 

~ 0.66 « 0.66 / t .0.66 x >0.66 

Po ” A \£) 

where p o 0 47 and p 0 °' 66 are the path radiances at 0.47 and 0.66 pm, respectively. The spectral 
dependence of the path radiance distinguishes between dust (dominated by coarse mode) and non-dust 
(dominated by fine mode). Kaufman, et al. (1997a) describes how the algorithm uses the ratio of path 
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radiances, p o 066 / p o 047 to make a 3 -branched decision whether the aerosol is pure dust, non-dust, or 
mixed. Assuming the continental model does not impact the resulting ratio. It is simply used to 
remove the molecular and surface contributions and to isolate the aerosol reflectance in both channels. 
The thresholds for the decision tree are: 

Po °' 66 / Po 047 < 0.72 THEN pure non-dust (3a) 

p 0 °' 66 / Po ° 47 > 0.9-0.0 1 (0 -1 50°) THEN pure dust (3b) 

0.72 < p 0 °' 66 / p o 047 < 0.9-0.01(0 -150°) THEN mixed (3c) 

for scattering angles 0= 150° to 168°. For scattering angles < 150°, 0 is simply, set to 150° for the 
boundary to collapse to 0.9 in that angle range. If the aerosol is mixed, then the fraction that the fine 
mode contributes to the total optical thickenss, T) = %d%tot is given by 


7 ] = 1 -- 


A 47 


-0.72 


0.90 - 0.0 1(0 - 1 50^) - 0.72 


(4) 


again, where 0 is set to 150° when 0 < 150°. Figure 2 shows four monthly mean values of the 
fraction of total aerosol optical thickness attributed to the 'non-dust' aerosol model. Red shades 
indicate that 'non-dust' dominates over the monthly mean. Purple shades indicate that 'pure dust' 
dominates. Blank areas in black are where no retrievals were made due to too bright a surface, 
monthly domination by clouds or snow. Note that these plots do not differentiate between high and low 
aerosol loading. Sensitivity to aerosol size decreases in very clean regions. 

Note that 'non-dust' is a misnomer because even though the non-dust aerosol models are 
dominated by their fine modes, each also contains a coarse mode as well. Three non-dust models are 
available and are described in Table 1, along with a description of the Continental model and the Dust 
model. The Urban/industrial model remains unchanged from Kaufman, et al., (1997b) and Remer and 
Kaufman (1998). The new Developing World - Moderate Absorption model is based on the 
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Biomass Burning model of Kaufman, et al. (1997b) and Remer et al. (1998), slightly modified to 
reflect the more recent study of Dubovik et al. (2002a). The Developing World - Heavy Absorption 
model uses the same size parameters as the other Developing World Model, but allows for the greater 
aerosol light absorption noted in Africa (Ichoku et al, 2003; Eck et al. 2003; Dubovik, et al. 2002a). 
Similar to the original conception of the algorithm, the current version uses season and geography to 
choose between the three non-dust models (Dubovik et al. 2002a). However, the boundaries have 
changed. Figure 3 gives the new distribution of the three non-dust models. 

If the aerosol is identified as dust by Equation 3, then the dust model of Table 1 is used. 
However, pure dust poses a problem with the dark target method. The assumption that p 2 ' 13 is 
transparent to aerosol and provides direct information from the surface does not hold when the aerosol 
is composed of large particles. Therefore, Equation (1) is not expected to hold in the pure dust case for 
very dark surfaces. However, over moderately bright surfaces, near the point of critical reflectance ( 
Kaufman, 1989) the surface contribution is negligible and the procedure can continue with minimal 
uncertainty introduced from the surface. Therefore, in the pure dust case retrievals are made only 
when p 213 falls between 0.15 and 0.25. 

In both the non-dust and pure dust cases, the estimated surface reflectances (p s ° 47 , p s 0 ' 66 ) and the 
measured mean top-of-atmosphere reflectances (p 047 , p 066 ) are used as input into the chosen model's 
Look-Up Table (LUT) to retrieve values for the aerosol optical thicknesses, fluxes and other 
parameters. A full second retrieval is made from the appropriate model's LUT. This second full 
retrieval differs from Kaufman et al. (1997b) that describes, instead, a correction based on the single 
scattering approximation to the preliminary Continental model retrieval described above. The final 
step in the process is to interpolate the values at 0.47 and 0.66 pm using an Angstrom law in order to 
report optical thickness and flux values at 0.55 pm. 
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The standard dark target retrieval path, described above, that uses specific aerosol models 
requires a minimum of 12 dark pixels in every 20 by 20 pixel nominal 10 km box. Path B, of Figure 1, 
represents an alternative if the surface is too bright to support 12 dark pixels in the standard manner. 
The upper limit of the p 213 value is permitted to increase as a function of the slant path. When the sun 
is at zenith and the satellite-view nadir, Path B collapses back to Path A requirements. However, as the 
photon path increases, more and more signal originates from the atmosphere, and the contribution from 
the surface reflectance becomes less and less important. This is especially true at the 0.47 pm 
channel, where atmospheric signal is highest and the surface usually darkest. In Path B, we only 
retrieve at 0.47 pm, for this reason. A final upper bound of p 2 13 = 0.40 is imposed. At least 12 pixels 
must again meet the Path B criteria, otherwise the procedure ends with no retrieval made, and fill 
values are placed in the output fields. Path B is considered to be less accurate than Path A, and the 
Quality Control is set to 0, ‘poor quality’. Because of the greater uncertainty over these brighter 
surfaces and because we retrieve in only one wavelength and cannot use the path radiance ratio to 
distinguish between dust and non-dust aerosol, only the Continental Model is used in the retrieval. The 
aerosol optical thickness and flux are derived from the LUT for 0.47 pm. These parameters are 
extrapolated to 0.55 pm and 0.66 pm using the spectral dependence of the Continental Model. 

Figure 4 shows a scene from the eastern part of southern Africa, where the surface reflectance 
is moderately bright. The top panel shows the MODIS retrieval of aerosol optical thickness at 0.55 pm 
following Path A of the algorithm, which insists on 12 dark target pixels. The bottom panel shows the 
modified version of the algorithm, which allows retrievals over brighter surfaces by following both 
Path A and Path B. From the image we see how the extension to brighter surfaces fills in holes without 
introducing suspicious artifacts. In this example, permitting Path B increases the number of retrievals 
over land from 7060 to 17,849. As we extend to brighter surfaces, we move away from the biomass 
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burning regions into cleaner regions. Thus the mean optical thickness of the granule decreases from 
0.20 to 0.15, but the standard deviation of the optical thickness remains constant at 0.15. For the 285 
granules collected over southern Africa during the SAFARI 2000 campaign, extending to brighter 
surfaces increases the number of land retrievals by 130%. 

2.2 The Ocean Algorithm 

The mechanics of the ocean algorithm are illustrated in Figure 5. Although the core inversion 
remains similar to the process described in Tanre, et al. (1997), the masking of clouds and sediments, 
the special handling of heavy dust including dust retrievals over glint, and revisions of the look-up 
table are new. As in the land algorithm, after the water vapor, ozone and carbon dioxide corrections 
are applied, the first step in die ocean algorithm is to organize the reflectance from the six wavelengths 
used in the procedure (p 0 55 , p 066 , p 0 ' 86 , p 124 , p 1 6 , p 213 ) into nominal 10km boxes of 20 by 20 pixels at 
500 m resolution. This requires degrading the resolution of the 250 m channels (p 0,66 and p 086 ). The 
ocean algorithm requires all 400 pixels in the box to be identified as ocean pixels by the MOD/MYD35 
mask. This helps to minimize problems introduced by shallow water near the coasts. If any land is 
encountered, the entire box is left for the land algorithm. 

Cloud and Sediment Masking: If all 400 pixels in the box are identified as water pixels, the 
algorithm then begins the arduous task of separating 'good' pixels from 'cloudy' pixels. The standard 
MOD35 cloud mask includes using the brightness in the visible channels to identify clouds. This 
procedure will mistake heavy aerosol as 'cloudy', and miss retrieving important aerosol events over 
ocean. On the other hand, relying on IR-tests alone permits low altitude, warm clouds to escape and be 
misidentified as 'clear', introducing cloud contamination in the aerosol products. Thus, our primary 
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cloud mask is based on the difference in spatial variability between aerosols and clouds (Martins, et al. 
2002). The algorithm marches through the 10 by 10 box, examining the standard deviation of p 0,55 in 
every group of 3 by 3 pixels. Any group of 9 pixels with standard deviation greater than 0.0025 is 
labeled as 'cloudy', and all 9 pixels in the group are discarded (Martins, et al. (2002)). The only 
exception to this rule is for heavy dust, which may at times be as spatially inhomogeneous as clouds. 
Heavy dust is identified by its absorption at 0.47 pm using the ratio (p 0 47 / p 0,66 ). This quantifies the 
difference that our eyes witness naturally. Dust absorbs at blue wavelengths and appears brown. 
Clouds are spectrally neutral and appear white to our eyes. If p 0 ' 47 / p 0 ' 66 < 0.75, then the central pixel of 
the group of 9 is identified as 'dust' and will be included in the retrieval even if it is inhomogeneous. 
This is a conservative threshold that requires very heavy dust in order to avoid clouds. Less restrictive 
thresholds would permit more dust retrievals, but might accidentally permit cloud contamination. 

The spatial variability test separates aerosol from most cloud types, but sometimes fails at the 
centers of large, thick clouds and also with cirrus, both of which can be spatially smooth. The centers 
of large, thick clouds are very bright in the visible, and so we identify these clouds when p° 47 > 0.40. 
This is an extremely high threshold that translates into an aerosol optical thickness greater than 5.0, but 
only for non-absorbing aerosol. Absorbing aerosol never reaches that high value of reflectance and 
will pass this cloud test unscathed. Some high values of non-absorbing aerosol may be discarded 
along with bright clouds, but this confusion is rare. Most heavy aerosol loading, with x > 5.0, absorbs 
somewhat at 0.47 pm and fails to reach the 0.40 threshold value, exhibited by very bright white clouds. 

Cirrus clouds are identified with a combination of infrared and near-infrared tests. Three 
infrared tests provided by the standard MODIS cloud mask, MOD35, are examined. These tests are IR 
cirrus test (byte 2, bit 4 ), 6.7 pm test (byte 2, bit 8) and Delta IR test (byte 3, bit 3) (Ackerman, et al. 
1998). If any one of the 3 indicates cloud, we label the pixel as 'cloudy'. The near-infrared cirrus test 
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is based on the reflectance in the 1.38 (am channel and the ratio p 138 / p 124 (Gao, et al. 2002). It is 
applied in the algorithm as a three step process. If p 138 / p 124 > 0.3, then the pixel is 'cloudy'. If 0. 1 0 
< p !38 / p 124 < 0.30 and p 138 > 0.03 and p 066 > l.Sp 0 ' 66 ^,^, then the pixels is also 'cloudy'. However, 
if 0.10 < p 138 / p 124 < 0.30 and 0.01 < p 138 < 0.03 and p 066 > l.Sp 0 66 ^^, then the situation is 
ambiguous. The algorithm labels the pixel as 'not cloudy' and will include the pixel in the retrieval 
process, but the quality of the retrieval (QC) is reduced to 0, 'poor quality'. This permits a retrieval at 
the orbital level (Level 2), but prohibits the retrieval from contributing to the long-term global aerosol 
statistics (Level 3). Only retrievals with QC > 0, contribute to the Level 3 Quality Weighted products. 
The products and product levels will be explained further in Section 3. Any situation in which the 
reflectance at 0.66 pm (p °’ 66 ) does not exceed 1.5 times 
the Rayleigh reflectance in that channel (p° 66 Rayleig h ) or the reflectance at 1.38 pm does not exceed 0.01, 
then the pixel is assumed to be 'not cloudy' with no ambiguity, unless the ratio (p 138 / p 1 24 ) exceeds 
0.3. 

The final mask applied to the data is the sediment mask, which identifies which ocean scenes 
are contaminated by river sediments (Li et al, 2002) and discards those pixels. The sediment mask 
takes advantage of the strong absorption by water at wavelengths longer than 1 pm. The resulting 
spectral reflectances over water with suspended sediments thus show elevated values in the visible, but 
not in the longer wavelengths. This creates an unique spectral signature quite different from clear 
ocean water and also different from airborne dust. 

The algorithm sorts the remaining pixels that have evaded all the cloud masks and the sediment 
mask according to their p 0 ' 86 value, discards the darkest and brightest 25%, and thereby leaves the 
middle 50 % of the data. The filter is used to eliminate residual cloud contamination, cloud shadows, 
or other unusual extreme conditions in the box. Because the ocean cloud mask and the ocean surface 
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are expected to be less problematic than their counterparts over land, the filter is less restrictive than 
the one used in the land retrieval. Of the 400 pixels in the original box, at least 10 must remain from 
the masking and filtering. Otherwise, no retrieval is attempted and fill values are given for the aerosol 
products in thatlO km box. If there are at least 10 good pixels in the 0.86pm channel and at least 30 
good pixels in the remaining 5 channels, the mean reflectance and standard deviation are calculated for 
the remaining 'good' pixels at the six pertinent wavelengths. 

Ocean Glint and Internal Consistency: The ocean algorithm was designed to retrieve only 
over dark ocean, away from glint. There is a special case when we retrieve over glint, and that is 
described below. The algorithm calculates the glint angle, which denotes the angle of reflection, 
compared with the specular reflection angle. The glint angle is defined as 
©g Im = co s _1 ((cosO s cos0 v )+(sinO s sin0vCOS(j))) (6) 

where 0 S , 0 V , and <() are the solar zenith, the satellite zenith and the relative azimuth angles (between the 
sun and satellite), respectively (Levy et al. 2002). Note that Fresnel reflection 
corresponds to 0 glint = 0. If 0 glint > 40°, we can avoid glint contamination and proceed with the 
retrieval. The algorithm performs several consistency checks of the spectral reflectances. Depending 
on the outcome of these consistency checks, the algorithm may either declare the reflectances to be 
beyond the range necessary for a successful inversion and exit the procedure, or continue onto the 
inversion after assigning quality flags (QC values) to each wavelength. 

Inversion Procedure: The inversion procedure is aptly described in Tanre et al. (1997) and 
Levy et al. (2002). Following Tanre et al. (1996), we know that the 6 reflectances measured from 
MODIS and used in the ocean retrieval (0.55 - 2.13 pm) contain three pieces of information about the 
aerosol. From this information we derive three parameters: the optical thickness at one wavelength 
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(x tot °' 55 ), the reflectance weighting parameter at one wavelength (p 0S5 ) and the effective radius, which is 
the ratio of the 3rd and 2nd moments of the aerosol size distribution. The inversion is based on a look- 
up table (LUT) that now consists of four fine modes and five coarse modes (Table 2, following Levy et 
al. (2002), which differs from the 11 possible modes listed in Tanre et al. (1997)). The LUT is 
constructed using the radiative transfer code of Ahmad and Fraser (1982). It consists of the top of 
atmosphere reflectances in six wavelengths calculated for a variety of geometries, a rough ocean 
surface with non-zero water-leaving radiance only at 0.55 pm (p s 055 = 0.005) and several values of 
T tot ° 55 for each single mode aerosol model of Table 2. Note that the LUT is defined in terms of a single 
wavelength of optical thickness. However, the parameters of each of the single mode models define a 
unique spectral dependence for that model, which can be applied to the retrieved value of T tot 0 55 to 
determine optical thickness at other wavelengths. Table 3 gives the spectral dependence of extinction, 
asymmetry parameter, and single scattering albedo for each of the modes of Table 2. 

The procedure requires both a fine mode and a coarse mode for each retrieval. The modes from 
the LUT are combined using T| as the weighting parameter, 

Pmr X (V 055 )^ (C' 55 )+(l-fr)PoW- 55 ) (7) 

The meaning of Eq. 7 is that the spectral reflectance measured from the satellite that 
corresponds to the LUT value Plot 1 " (Lot 0 ' 55 ) for the determined values of p and i; ot 0 ' 55 is a weighted 
average of the reflectance values for an atmosphere with a pure fine mode 'f and optical thickness 
t; tot 0 55 and the reflectance of an atmosphere with a pure coarse mode 'c' also with the same T tot ° 55 . The 
weighting factors are T) and 1- p. In Appendix A, we show that p = p 0 55 , the fraction of total optical 
thickness at 0.55 pm contributed by the fine mode. 
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For each of the twenty combinations of one fine mode and one coarse mode, the inversion finds 
the pair of T^ 0 ' 55 and p 0 ' 55 that minim izes the error (e) defined as 

2 

( 8 ) 


where N x is the sum of good pixels at wavelength X, p m x is the measured MODIS reflectance at 
wavelength X, and Pm/ is calculated from the combination of modes in the look-up table and 
defined by Equation (7). The 0.01 is to prevent a division by zero for the longer wavelengths under 
clean conditions (Tanre et al. 1997). The inversion requires P L ut°' 87 to exactly fit the MODIS 
observations at that wavelength and then finds the best fits to the other five wavelengths via Equation 
(8). The 0.87 pm channel was chosen to be the primary wavelength because it is expected to be less 
affected by variability in water leaving radiances than the shorter wavelengths, yet still exhibit a strong 
aerosol signal, even for aerosols dominated by the fine mode. By emphasizing accuracy in this channel 
variability in chlorophyll will have negligible effect on the optical thickness retrieval and minimal 
effect on p° 55 . 

The twenty solutions are then sorted according to values of e. The best solution is the 
combination of modes with accompanying x tot 0 55 and p 0S5 that mini mizes £. The solution may not be 
unique. The average solution is the average of all solutions with £< 3% or if no solution has £ < 3%, 
then the average of the 3 best solutions. Once the solutions are found, then the chosen combination of 
modes is the de facto derived aerosol model and a variety of parameters can be inferred from the 
chosen size distribution including spectral optical thickness, effective radius, spectral flux, mass 
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concentration, etc. 
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Final Checking. Before the final results are output, additional consistency checks are 
employed. In general, if the retrieved optical thickness at 0.55 pm is greater than -0.01 and less than 

5, then the results are output. Negative optical depths are given lower quality flags. There are 

I 

exceptions and further checking for heavy dust retrievals made over the glint. Quality flags may be 
adjusted during this final checking phase. 

Special case: Heavy dust over glint. If 0 glint < 40° then we check for heavy dust in the glint. We 
use a similar technique as before during the masking operations when we noted that heavy dust has a 
distinctive spectral signature because of light absorption at blue wavelengths. In the situation of 
identifying heavy dust over glint we designate all values of p 047 / p° 66 < 0.95 to be heavy dust. If 
heavy dust is identified in the glint, the algorithm continues with the retrieval, although it sets QC=0. 
This permits the retrieval, but prohibits the values from being included in the Quality Weighted Level 
3 statistics. If heavy dust is not identified in the glint, then the algorithm writes fill values to the 
aerosol product arrays and exits the procedure. 

3.0 The Aerosol Products 

| 

Examples of the three main aerosol products are shown in Figure 6. These three main products are 

I 

the land and ocean aerosol optical thickness at 0.55 pm (t 0 55 ), the fraction of the optical thickness 
contributed by the fine mode (q 0 55 ) and the reflected flux at the top of the atmosphere at 0.55 pm. The 
flux products are hemispherical irradiances at a particular wavelength. Over land, the flux is defined 
for zero surface reflectance and computed consistently with the same aerosol parameters used in the 
optical thickness retrieval. Over ocean, the flux is defined for the same rough ocean surface model and 
the same aerosol parameters derived in the optical thickness retrieval. Errors introduced by inaccurate 


i 
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assumptions of the aerosol model during the retrieval when measured radiances at top of atmosphere 
are converted to optical thickness are subsequently canceled when those same aerosol model 
assumptions are used with the optical thickness to calculate top of atmosphere fluxes. Errors 
introduced by assumptions of surface reflectance will remain. Over a monthly average, the consistent 
calculations produce a highly accurate measure of the flux, more accurate than the optical thickness 
itself. Figure 6 shows orbital-based Level 2 products, which are cut from 5 minute segments of a 
satellite orbit, called granules. The Level 2 products, designated as MOD04 files, contain 64 separate 
products, all connected to the aerosol retrieval. The Level 2 files are produced every day and represent 
the first level of MODIS aerosol retrieval. In addition, statistics based on the Level 2 aerosol 
retrievals can also be found in Level 3 files, designated as MOD08 files. These Level 3 files contain 
parameters produced from the entire MODIS Atmospheres team and include such parameters as water 
vapor and cloud characteristics along with the aerosol information. The Level 3 data are averaged to a 
1 -degree latitude/longitude grid and are produced every day (MOD08_D3), averaged every 8 days 
(MOD08_E3) or averaged on a monthly basis (MOD08_M3). They include both statistics calculated 
equally from all the data, and also statistics weighted by the quality of each individual retrieval. 
Quality weights of 0 will prevent poor retrievals from affecting the calculated statistics of the Quality 
Weighted quantities. Further information about the Level 3 products can be found in (King, et al. 
2003)) and at http://atmos-modis.gsfc.nasa.gov . 

All together there are 64 aerosol products at Level 2: 9 products describe geometry and 
location, 3 products are joint land and ocean products, 23 are land-only products and 29 are ocean-only 
products. Tables 4-7 lists all 64 products. The three joint land and ocean products are simple 2- 
dimensional arrays of one wavelength (Figure 6). The land only and ocean only products contain an 
additional dimension. In many cases this additional dimension is wavelength. Tables 4-7 list the 
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wavelengths for each product, where applicable. The additional dimension in the ocean-only products 
can designate either the 'best' solution or the 'average' solution from the ocean retrieval as described 
above in Section 2.2 (Inversion Procedure). Both solutions are reported for some parameters, although 
they are often identical. 

Tables 4-7 also list whether the product is 'validated', 'not yet validated', 'derived', 'experimental', or 
'diagnostic'. A 'validated' product indicates that substantial comparison was made to ground-based 
data, and that the retrieval is well-characterized so that error bars can be defined and comfortably 
applied to the retrieval product (Ichoku et al.,2002; Chu et al., 2002; Remer et al. 2002). Examples of 
validation will be shown in Section 4. ‘Not yet validated’ indicates that the retrieved parameter has not 
yet been well-characterized, but that data is being collected and analysis is underway. ‘Derived’ is a 
parameter that follows from the retrieval’s choice of aerosol model and the magnitude of the retrieved 
optical thickness. Definitions of some of the derived parameters are given in Appendix B. A ‘derived’ 
parameter is not directly retrieved and there are no expectations of ever validating a ‘derived’ 
parameter with independent data. ‘Experimental’ is a scientific product that may have future 
applications, but as of now is too innovative to be well-characterized. ‘Diagnostic’ refers to output that 
is either an auxiliary or intermediate parameter. ‘Diagnostic’ parameters are meant to aid in 
understanding the final product, but will never themselves become ‘validated.’ 

Recommendations for choosing particular products are given in Appendix C. 

4.0 Validation of Aerosol Products 

Our primary means of validation is comparison with AERONET ground radiometers (Holben, 
et al. 1998). The AERONET instruments measure spectral aerosol optical thickness, to within 
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— 0.01 in the MODIS wavelengths (Eck et al., 1999) and can derive ambient, total atmospheric column 
aerosol effective radius, r_eff, whenever conditions are favorable (Dubovik, et al. 2000). The 
methodology of comparing temporally varying AERONET data with spatially varying MODIS data is 
described in (Ichoku et al., 2002). In the following validation, we use AERONET Level 1.5 data, 
which are cloud screened but not quality assured, primarily because final calibration is not applied 
(Smirnov, et al. 2000). The procedure that co-locates MODIS and AERONET data is applied during a 
very short window of opportunity (~24 hours) after the operational MODIS aerosol product has been 
processed and before that data is transferred to the archive. The data is purged from the operational 
facility after transfer. It would be a monumental task to retrieve the global data back from the archive 
to re-link with the AERONET Level 2.0 quality assured data. Figure 7 shows the distribution of the 
132 AERONET stations used in the comparisons to be described below. Although North America and 
Europe provide the most stations in the data base, all continents (except Antarctica), all oceans and all 
aerosol types are represented. 

Validation is an on-going effort. Not only do aerosol conditions vary in location and time, 
requiring a continued effort to validate the algorithms under various conditions, but the algorithms 
themselves evolve. The algorithms development and history, starting from the most recent version and 
going backwards in time can be found at the MODIS atmospheres web site (http://modis- 
atmos.gsfc.nasa.gov/MOD04_L2/history.html). Much of the algorithm modifications concern 
experimental products , improved cloud /snow/water masking, internal bookkeeping or minor 
adjustments meant to improve long-term statistics. However, the last adjustment over land to Version 
4.1.0 modifies the land look-up tables to increase aerosol absorption as needed in certain regions 
(Ichoku et al. 2003). This is a significant change, which will not be reflected in the validation plots 


described below. 
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A preliminary validation of the aerosol products was made of the data collected in the first 
months of operation. The results are reported in Ichoku et al. (2002), Chu et al. (2002) and Remer et 
al. (2002). The preliminary validation compared 2 to 3 months of MODIS aerosol optical thickness 
and effective radius retrievals to the same parameters observed (optical thickness) or derived (effective 
radius) from AERONET radiometers. The preliminary validation from the limited data set showed 
good agreement between the MODIS-derived parameters and the AERONET parameters. 

4.1 Validation of Aerosol Optical Thickness 

Validation Over Land. Currently there are two years of MODIS aerosol products co-located with 
AERONET retrievals (August 1, 2000 to August 1, 2002). Figure 8 shows plots of 5906 co-located 
points over land at wavelengths 0.47 pm, 0.55 pm and 0.66 pm. There are no AERONET 
measurements at the MODIS wavelengths of 0.47 pm and 0.55 pm; therefore, the AERONET values 
in the plots of Figure 8 have been interpolated from the values at 0.47 pm and 0.87 pm. The 0.50 pm 
AERONET channel is not used for interpolation because not all AERONET stations have that channel 
and the procedure that matches MODIS and AERONET data must be uniform and automatic. 
Although AERONET does make measurements at 0.675 pm, the values at this wavelength have also 
been interpolated from 0.44 pm and 0.87 pm, due to occasional calibration drift at this channel in the 
AERONET instruments. These calibration issues due to gradual filter degradation are identified and 
corrected in the Quality Assured Level 2 AERONET data, but unfortunately the co-location with 
MODIS is done in real-processing time and cannot wait for the post-deployment AERONET 
calibration corrections. Therefore, the 0.675 pm is not used and instead the information is transferred 
from the more reliable 0.44 pm and 0.87 pm AERONET channels. Also, the MODIS value at 0.55 
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um is not a direct retrieval, but an interpolation from the 0.47 (am and 0.66 pm retrievals; thus, the 
plot at 0.55 pm in Figure 8 is a comparison of two interpolated values. 

Figure 8 represents the scatter plot between MODIS retrievals and AERONET observations, 
co-located in space and time. The data were sorted according to AERONET aerosol optical thickness. 
Then an average was calculated for every 300 points and plotted. At higher optical thickness where 
the data become sparser, fewer points are used in the average, as indicated. The standard deviation in 
each bin is shown by error bars. The regression equation and correlation given at the top of each plot 
were calculated from the full scatter plots, before binning. The solid black line is the 1:1 line, and the 
dashed lines denote the expected uncertainty calculated from pre-launch analysis. These dashed lines 
should encompass one-standard deviation (66%) of the aerosol retrievals. The pre-launch expected 
uncertainty over land is Ax = ± 0.05 ± 0.15t (Chu et al., 1998; King, et al. 1999). The regression 
equations in Figure 8 show all wavelengths having an offset greater than the expected offset of 0.05 at 
low optical thickness, while the plots show a positive bias at low optical thickness suggesting that there 
is an issue with instrument calibration or more likely, that surface reflectance may be improperly 
represented in a systematic way at certain locations and seasons. The regression equations also show 
all wavelengths with slopes less than one. Ichoku et al. 2002 demonstrate that this underprediction of 
aerosol optical thickness at higher aerosol loadings can be attributed to underpredicting the amount of 
light absorption by the aerosol in certain regions of the world. This is the reason why the Strong 
Absorption Model (Table 1) was introduced to the algorithm in the Version 4 delivery. We expect the 
under prediction of optical thickness at high aerosol loading to be less of an issue after reprocessing 
with the updated algorithm, but the issue of the offset at low aerosol loading will remain. 

Even with the deviations described above, the results of Figure 8 indicate that the algorithm is 
retrieving aerosol optical thickness over land to roughly within the expected accuracy. On a global 
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basis, 61%, 68% and 71% of the retrievals at the 0.47 urn, 0.55 pm and 0.66 (am, respectively, fall 
within expected error. From a global perspective, only the blue channel is falling outside of the error 
bars slightly more often than the pre-launch expectations of 66%. The average x at 0.55 pm is 0.1 8 for 
the land global data base at defined AERONET stations and the percent error between MODIS 
retrievals and AERONET observations at 0.55 pm is 41%, showing a positive bias in which MODIS 
overestimates x. The overestimate corroborates the positive offsets seen at low to moderate values of 
optical thickness in Figure 8. Table 8 shows the percent of retrievals falling within the expected error 
lines for all the data, and also grouped by specific region. In some regions, the retrievals are more 
consistently accurate than in other regions. Specifically the North American continent, especially 
Alaska/Canada is proving to be difficult. This region also exhibits the lowest x, which contributes to 
the high relative error. The absolute error in Alaska/Canada is comparable to other regions. 

Validation Over Ocean. Similarly to Figure 8 for the land retrievals, Figure 9 represents the scatter 
plot of 2052 MODIS retrievals over ocean co-located with an AERONET station either on the coast or 
on an island. Note the dashed lines, denoting expected uncertainty, are narrower than those over land. 
The MODIS over-ocean algorithm is expected to be more accurate than the over-land algorithm (Ax = 
± 0.03 ± 0.05x ), (Tanre et al., 1999; King et al., 1999). AERONET values at 0.55 pm and 0.66 pm, 
are interpolated as in Figure 8. No MODIS values are interpolated for these plots, and the 0.87 pm 
plot is the only one showing a directly retrieved MODIS value plotted against a directly measured 
AERONET value, with no interpolation for either quantity. Unlike the land validation of Figure 8, the 
ocean algorithm has virtually no offset and little bias, except for a possible absorbing aerosol bias at 
high optical thickness suggested in the 0.87 pm plot. The linear regression line follows the 1:1 line 
closely through out the optical thickness range where most retrievals occur. 
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Table 8 demonstrates the overall higher accuracy of the ocean retrieval when compared to land 
in that the percent error is consistently smaller over ocean than over land. Globally, 62%, 66% and 
70% of all retrievals over ocean at 0.55 pm, 0.66 pm and 0.87 pm, respectively, are falling within the 
narrowly defined expected uncertainty. Only the 0.55 pm channel is falling outside of the error bars 
slightly more often than the pre-launch expectations of 66%. The average T at 0.55 pm is 0.18 for the 
ocean global data base at defined AERONET stations, the same as for land. Because the land and 
ocean data bases include many of the same stations, this is not surprising. The percent error between 
MODIS ocean retrievals and AERONET observations at 0.55 pm is only 1%, showing the same 
absence of bias as exhibited in Figure 9. 

Regionally Table 8 shows that the AsianPacific_ocean region and the Saharan_ocean region 
fall outside the expected uncertainty lines more often than other regions, although for the Sahara it is 
random scatter with no preference as to over or under prediciting. The aerosols in these regions can 
have a strong dust component. (Levy et al., 2003) demonstrate that the ocean-algorithm does not 
perform well in a dust-laden aerosol, attributing the problem to poor assumptions for the dust-aerosol 
phase functions (ie. Nonsphericity). Empirical non-spherical phase functions have been derived and 
will be implemented into the next version of the MODIS ocean algorithm. We expect improvements 
for retrievals in dusty regions to follow from this update. 

Discussion of Validation. Comparison of MODIS retrievals with highly accurate ground-based 
radiometer data validates the basic retrieval, but does not necessarily validate the product for use in 
long-term climate studies. Figures 8 and 9 cannot validate the MODIS cloud clearing algorithms that 
play a significant role in the quality of the retrieved aerosol products. Figures 8 and 9 represent the co- 
located points only for the events in which both the MODIS and the AERONET cloud masking 
algorithms indicated that no clouds were present. It is possible that MODIS might retrieve, while 
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AERONET identified clouds and did not. Those points would not show up on the scatter plots. If 
those cases were numerous, long-term MODIS aerosol statistics could be cloud-contaminated despite 
the good agreement in figures 8 and 9. 

Figure 10 presents an alternative method of validation that checks the long-term statistics for . 
possible cloud contamination. Here we plot monthly mean values in a 3 degree latitude by 3 degree 
longitude box centered on the AERONET station. The data is for the year 2001 at eight selected 
stations, 4 with land components and 4 with only ocean retrievals. Within the 4 land sites, three of 
them are near enough to the coast to contain sufficient ocean retrievals within the 3 degree by 3 degree 
box and can be used for testing both land and ocean retrievals. The thick lines, in each plot represent 
the MODIS aerosol optical thickness at 0.55 pm for land (red) and ocean (blue), and the AERONET 
value (black). The MODIS monthly mean values were calculated from archived MODIS Level 3 data 
(MOD08) on a 1 degree resolution. Thus, the difficulties with matching MODIS with AERONET in 
near-real time do not exist, and the monthly AERONET values are calculated from AERONET Level 
2.0 data. The data from MODIS and AERONET need not be simultaneous. Despite the non- 
simultaneous nature of the measurements, both MODIS and AERONET are expected to represent 
monthly mean aerosol optical thickness for each region. In this comparison, MODIS does not benefit 
from AERONET's cloud clearing algorithm. If MODIS retrievals were systematically cloud 
contaminated, we would expect the MODIS monthly mean values to be systematically higher than 
AERONET's. This is not the case in any of the 8 regions examined. In most cases, MODIS and 
AERONET exhibit very similar annual cycles, often with very similar magnitudes of optical thickness. 
When magnitudes systematically differ, as in Cuiaba-Miranda of Figure 10, it is AERONET that 
systematically exceeds MODIS. The spatial variability across the 3-degree box can explain some of 
the differences between instruments. Still, the agreement in optical thickness is striking and suggests 
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that MODIS monthly mean optical thickness values are not significantly cloud contaminated and can 
be used with confidence in developing a global aerosol climatology and estimating aerosol forcing. 

4.2 Validation of Aerosol Size Parameters 

The validation of retrieved size parameters is not as straightforward as validation of optical 
thickness. Our primary means of validation is to compare with derivations of the same parameter from 
inversions of AERONET observed sky radiance ((Dubovik et al., 2000)). Sky radiance measurements 
are taken less often than direct sun measurements in the AERONET protocol. Furthermore, sky 
radiance data must be sufficiently homogenous and the inversion must make a good fit to the measured 
radiances in order for the retrieval to be used. As a result there are fewer simultaneous data to be 
plotted in a scatter plot. Because of this we rely primarily on comparisons of monthly means, which 
test the applicability of the long-term statistics. 

Figure 10 shows comparisons of monthly mean MODIS and AERONET derived r) 0 55 (the ratio 
of fine mode to total optical thickness). These are depicted as the thin lines in the figures, with blue a 
MODIS ocean derivation, red a MODIS land derivation and black from AERONET. The AERONET 
values are calculated from standard inversions of AERONET observed sky radiance (Dubovik et al., 
2000). The MODIS size parameters over land are not expected to be as accurate as the parameters 
over ocean, and therefore we focus our discussion on the ocean-derivations shown by the thin blue 
lines (Tanre et al. 1996). For some sites such as GSFC, Anmyon, and Male, MODIS-ocean and 
AERONET agree to within 20% for much of the year. For Bermuda, Midway _Island and Lanai, the 
agreement is sustained for the first 6 months of the year until MODIS-size parameter jumps to a much 
lower value. These latter ocean stations exhibit very low optical thickness. The derivation of size 
parameters at low optical thickness, when aerosol signal is small, will be very sensitive to instrument 
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calibration. In June 2001, the MODIS instrument suffered an anomaly and the data processing was 
switched from the B-side to the A-side electronics. This switch created a small aberration in the 
calibration that did not affect the more robust optical thickness retrievals, but did affect the size 
parameter derivations when optical thickness was low. 

Comparison of MODIS-derived particle effective radius and AERONET-derived particle 
effective radius is shown in Figure 11. This plot matches MODIS retrievals to daily averages of Level 
2 Quality Assured AERONET sky radiance inversions. The MODIS effective radius parameter is for 
the total bimodal size distribution in which each of the two modes is weighted by T| (Appendix B). 
There were 492 MODIS effective radius retrievals co-located with a daily average AERONET 
retrieval. However only the 271 match-ups with optical thickness greater than 0.15 are plotted in 
Figure 12. At low optical thickness there is greater susceptibility to all algorithmic and sensor 
uncertainties including small calibration errors (discussed above) and retrieval errors for both 
instruments (Ignatov et al. 1998, Remer et al., 2002). These errors make little difference in optical 
thickness retrieval but create large errors when size parameters are calculated. Figure 1 1 shows that 
62% of the points fall within the ±0.10 pm errors, which is a smaller percentage than what Remer et al. 
(2002) reported for a more limited data set. Note that the MODIS and AERONET retrievals both 
assume spherical particles in deriving size distribution. This assumption causes both MODIS and 
AERONET to underpredict particle effective radius when non-spherical dust is present (Dubovik et al. 
2002b). Thus, the agreement in Figure 1 1 may be better at some sites if AERONET retrievals had 
assumed spheroids instead of spheres (Dubovik et al. 2002b). The MODIS algorithm is being 
modified to include the option of empirical phase functions that do not require any assumption of 
particle shape. 
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5.0 Results 

Figure 12 illustrates the MODIS aerosol retrievals at the global scale. The images are constructed from 
the aerosol optical thickness and size parameter products derived from observed MODIS radiances. 
Red indicates aerosol dominated by small particles (less than 0.5 pm) and greenish tints indicate 
aerosol with a higher proportion of large particles (greater than 0.5 pm). We can see that aerosol from 
natural sources, such as sea salt and desert dust, contain larger particles than aerosols emanating from 
human-produced combustion sources such as agricultural and deforestation burning or urban/industrial 
pollution. Thus, we see that aerosol size, with the exceptions of lightning-initiated forest fire smoke 
and ocean DMS production easily separates aerosol into natural and man-made components. Thus we 
see that MODIS's ability to separate aerosols by size can be used as a proxy for separating 
anthropogenic aerosol from natural sources, and increases the accuracy of estimating human-induced 
aerosol forcing (Kaufman, et al. 2002). 

6.0 Conclusions 

Characterizing the global aerosol system is essential to understanding the earth's climate system 
and estimating potential global climate change. The MODIS instrument flying aboard NASA;s Terra 
and Aqua satellites provides a look at the aerosol system over both land and ocean on a daily basis. 
The derivation of aerosol products from the MODIS -measured radiances relies on the broad spectrum 
MODIS measures, ranging from the visible into the mid-infrared, and the 500 m spatial resolution, 
which allows for better cloud identification and clearing than was possible with previous instruments. 
The mature MODIS algorithm includes aerosol optical thickness at several wavelengths, information 
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on particle size, and aerosol reflected flux at the top of the atmosphere, which is expected to be more 
accurate than the optical thickness retrievals. An extensive validation effort that co-located over 8000 
MODIS retrievals with AERONET measurements of optical thickness show that globally, the MODIS 
products are accurate to within pre-launch expectations. In particular, the retrieval of aerosol over 
oceans consistently shows remarkably good agreement through the range of optical thickness where 
most observations occur. Regional analysis shows specific issues for certain locations. Comparison of 
MODIS and AERONET monthly means at eight specific locations scattered globally demonstrates that 
the MODIS retrievals are not affected by cloud contamination at those sites, and that MODIS long- 
term statistics agree well with AERONET and can be used to determine global aerosol climatology. 
MODIS-derived aerosol size parameters are in general agreement with the same quantities derived by 
AERONET instruments on the ground. For moderate optical thickness, one-standard deviation of 
MODIS effective radius retrievals falls within ±0.11 pm of AERONET measurements. However, at 
low aerosol optical thickness ( x < 0.15) the MODIS size retrievals are susceptible to small aberrations 
in the calibration and other factors which introduce greater uncertainty. In addition, dust, with its 
nonspherical shapes, introduce uncertainty in both the optical thickness and size parameter retrievals. 
This latter issue will be addressed with the incorporation of non-spherical phase functions into the next 
version of the algorithms. In the meantime, the MODIS aerosol products are sufficiently accurate for 
a variety of applications including improved estimates of observationally based aerosol optical 


thickness. 



31 


Appendix A: Fraction of fine mode optical thickness. 

One of the most important products produced by the algorithm is the ratio of fine mode optical 
thickness to the total optical thickness, or simply the fraction of fine mode. Here we show that this 
fraction at 0.55 pm is the same parameter as, r), the reflectance weighting parameter. We start with 
equation (7) 

Plot W 5 ' 55 )^ Pf X ( T tot°' 55 ) + ( 1_T l)Pc X (‘ T tot 0:5S ) (Al) 

where p f and p c are the fine and coarse mode atmospheric reflectances for the same optical 

X 

thickness as the total spectral reflectance, p LUT , respectively, and r\ is the 

X- 

reflectance weighting parameter. Note that p LUT * p f + p c and that the total and component 

reflectances all are functions of the total optical thickness (x tot °' 55 ) not the component optical 

thicknesses (x f 0 ' 55 and x c 0 ' 55 ). All optical thicknesses are defined at 0.55 pm. This is by definition in 
constructing the look-up tables. 

Define the total optical thickness (x tot 0 ' 55 ) equal to the sum of the fine (x f 0 ' 55 ) and coarse (x c 0 ' 55 ) 
components. Using the single scattering approximation, 

Pf X = Cx tot °' 55 P/, Pc X = Cx tot °' 55 P c X (A2) 

and 
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p LUT X =C(T f °- 55 P/+ X C °' 55 P C X ) (A3) 

where C is a constant depending on geometry and Pf X and P c x are the fine mode and coarse mode 

phase functions calculated for the look-up tables, respectively. There is no P tot x because the phase 

functions in the look-up table are calculated for the collection of individual fine and coarse modes, not 
for any ‘total’ aerosol size distribution. Solving for p in Equation A1 gives, 

= (Plut^ — PcW — Pc ) (A4) 

Substituting, A 2 and A3 into A4 gives 

p = (C(T f 0 - 55 P f X + T c 0 - 55 P C X ) - Cl tot °' 5S P c 1 )/ 

(Cx tot °- 55 P f X -Cx tot 0 S 5 P c X ) (A5) 

Dropping the constant ‘C’ and using the definition of x tot 0 ' 55 = x /' 55 + x c 0 ' 55 gives 

p = (x f 0 - 55 P f x + x c °' 5s P c x ) - x f 0 - 55 P c x - x c 0 - 55 P c x )/ 

[Xtot 0 ' 55 (Pf X -P c X )] (A6) 

T1 = Xf 0 ’ 55 (Pf X - P C X )/[T tot °- 55 (P f X - P C X )] (A7) 


= 0.55/ ( 

n Xf ^tot 


0.5S 


(A8) 
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Thus the reflectance weighting factor, T|, is also the ratio between fine mode and total optical thickness 
at 0.55 pm, as defined within the parameters of the inversion. 
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Appendix B: Definitions of Derived Parameters 

The following give the formulas for derivation of the derived parameters. In these formulas n(r) is the 
size distribution with r denoting radius. r g is the geometric mean radius. N 0 is number of particles per 
cross section of the atmospheric column (the amplitude of the number size distribution of a lognormal) 
that can be converted from V 0 of the volume size distributions for each lognormal mode using 


3aV 0 V27t _3 , 9 2\ 

N 0 = 4^ r g e*p(-- 0 ) 

(3s is the scattering coefficient, specific to each model mode, p is the density of the particle assumed to 
be 1 g/cm . erfQ is the error function, c is In a g where a g is the geometric mean standard deviation of 


the lognormal distribution, p is cos 9, and 0 is the scattering angle. P() is the phase function, 
x is the optical thickness and unless designated specifically for wavelength or large or small mode, 
represents the total optical thickness at 0.55 pm. Extcoeff is the extinction coefficient and unless 
designated specifically for wavelength is understood to be 0.55 pm. 


Cloud condensation nuclei in units of cm' 2 , 
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v j J 


CCN = f°° . n(r)dr = 0.5x10' 

• l r o =0.03pm v ’ 

Assymetry factor, 

n = OS? hQlQMi = 0.5f* cos(d)/%(9)sm9Jd 


Backscattering ratio, 
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P\ = IT f'.-PftOdti = J-J^ 9 P( 9 )sm 9 de 

27C J — 1(J, 271 J{ J 


Number of particles in each mode (1 particle/cm 3 ), x and extcoeff defined at 0.55 pm, 
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Mass concentration (pg/crrT), 
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i Appendix C: Recommendations For Using Products: 

There are many choices for aerosol optical thickness. The products In Tables 4-7 labeled as 
‘validated’, ‘not yet validated’, or ‘derived’ are recommended. 'Corrected_Optical_Depth_Land' is the 
recommended spectral product over land, and 'Effective_Optical_Depth_Average_Ocean' is the 
recommended spectral product over ocean. Products such as ‘Continental_Optical_Depth_Land’ are 
intermediate parameters and should be used only as a diagnostic. 

The word 'small' in the product name indicates fine mode so that 'Optical_Depth_SmaH' is the fine 
mode optical thickness, T f , from Appendix A, 'Optical_Depth_Ratio_Small' is the fine mode ratio, T|\ 
Likewise, the word ’large' indicates coarse mode. The word 'Average' indicates the solution averaged 
from all retrieval solutions with fitting error less than 3% or the average of the 3 best solutions if all £ 

■ 

> 3 %. The word 'best' indicates the single solution with the least error, e, no matter how large. See 
equation 7. The recommendation is to use those products labeled as 'average'. 

I 

The Angstrom Exponent over land is defined for wavelengths 0.47 and 0.66 pm. There are two 
Angstrom Exponents for the ocean parameters, one defined using wavelengths 0.55 and 0.87 pm and 
the other using 0.87 pm and 2.13 pm. 

‘Aerosol_Type’ under the land products is a function mostly of geography and season and should 
not be considered a retrieved quantity. The ‘Cloud_Fraction’ listed in the tables is not a true cloud 
fraction, but instead an indication of the fraction of pixels not used in the retrieval due to a 
combination of clouds, surface issues or internal inconsistencies. Likewise, ‘Mean_Reflectance’ is the 
mean reflectance only of those pixels that survive the masking and elimination procedures and are 
actually used in the retrievals. Solution_Index tells which fine and coarse aerosol models were chosen 
in the retrieval. Least_Squares_Error reports the fitting error of the inversion, e, from Equation 7. 
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The ‘Quality_Assurance’ parameters are 5-byte codes that hold information concerning the 
retrievals and the overall quality. Details of the ‘QualityAssurance’ code are given by the MODIS 
Atmosphere’s Quality Assurance Plan that can be found at http://modis-atmos.gsfc.nasa.gov . The top 
bit of the Quality Assurance code is the most important. If this top bit equals 0, then the retrieval at 
that location is experimental. It may be set to 0 because the surface was bright and Path B of Figure 2 
was followed, or it may be a special heavy dust retrieval over ocean. Mostly these experimental 
retrievals provide excellent results, but sometimes they do not, and the user should be cautious when 
using data with Quality of 0. 
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Table Captions 

Table 1. Size Distribution Parameters, and Single Scattering Albedo Used in the MODIS Look Up Table 
for the Land Algorithm. 


Table 2. Refractive Indices, Median, Standard Deviation and Effective Radius for the aerosol models used 
in the MODIS Lookup Table for the ocean algorithm. Models 1-4 are fine modes and Models 5-9 are 
coarse modes. 


Table 3. Values of Assymetry parameter and Single Scattering Albedo for the 9 ocean models of Table 2. 
Table 4. Contents of MODIS Aerosol Level 2 hdf file (MOD04): Time and geometric information 
Table 5. Contents of MODIS Aerosol Level 2 hdf file (MOD04): Global land and ocean products, at 0.55 prr 

Table 6. Contents of MODIS Aerosol Level 2 hdf file (MOD04): Land products 

Table 7. Contents of MODIS Aerosol Level 2 hdf file (MOD04): Ocean products 

Table 8. Number of retrievals (N), percentage of retrievals (%) falling within expected uncertainty for each ( 

wavelength, average optical thickness from AERONET (t 550 ) and percent difference between MODIS and 

AERONET at 0.55 pm (Diff) for the global data set and for each region, land and ocean separately. 
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Figure Captions 


Figure 1. Flowchart illustrating the derivation of aerosol over land. 

Figure 2. Monthly mean plots of fraction of total aerosol optical thickness attributed to 'non-dust' or fine 
mode aerosol over land. Fraction 1.0 indicates all fine mode. Fraction of 0.0 indicates all coarse mode. 

Figure 3. Distribution of the non-dust models used in the derivation of aerosol over land. Single 
scattering albedo values given in parentheses. Uncolored regions use the urban/industrial aerosol model. 
Solid black regions use the moderate absorption aerosol model all year. Large checkerboard pattern in east 
Asia denotes a region of strong absorption aerosol model all year. Africa is divided by region and season. 
North of the equator, during the burning season (November-May) the strong absorption aerosol model is 
used, while the moderate absorption model is used the remainder of the year. South of the equator, the 
burning season shifts to June-October when the strong absorption model is used, while in the remainder of 
the year the algorithm uses the urban/industrial model. 

Figure 4. MODIS-derived aerosol optical thickness at 0.55 pm for an image of the east coast of southern 
Africa. Top panel uses the traditional dark target method described by Path A in Figure 2. The bottom 
panel shows the results after extending the retrieval to brighter surfaces as described by Path B of Figure 
2. By extending to brighter surfaces the number of retrievals over land in this image increases from 7060 
to 17,849. 


Figure 5. Flowchart illustrating the derivation of aerosol over ocean. 
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Figure 6. Examples of MODIS aerosol products at the Level 2 stage (MOD04). The data represent a 5- 
minute granule collected on July 7, 2002 from 1835-1840 UTC. The upper left panel is a true color image 
created from Level lb reflectances. The upper right panel is aerosol optical thickness at 550 nm, The lower 
left panel is reflected flux at 550 nm. The lower right panel is fraction of fine mode. 

Figure 7. The distribution of the 132 AERONET stations used to validate MODIS land and ocean aerosol 
retrieval algorithms. 

Figure 8. MODIS aerosol optical thickness retrievals over land at 470 nm (blue), 550 nm (green ) and 660 
nm (red) as a function of AERONET observations co-located in space and time. The data were sorted 
according to AERONET aerosol optical thickness and averaged for every 300 points. At higher optical 
thickness where the data become sparser, fewer points are used in the average, as indicated. The standard 
deviation in each bin is shown by error bars. The regression equations given at the top of each plot were 
calculated from the full scatter plots before binning. The solid black line is the 1 : 1 line and the dashed lines 
denote the expected uncertainty calculated from pre-launch analysis. 

Figure 9. MODIS aerosol optical thickness retrievals over ocean at 550 nm (green ) and 660 nm (red) and 
870 nm (black) as a function of AERONET observations co-located in space and time. The data were 
sorted according to AERONET aerosol optical thickness and averaged for every 100 points. At higher 
optical thickness where the data become sparser, fewer points are used in the average, as indicated. The 
standard deviation in each bin is shown by error bars. The regression equations given at the top of each 
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plot were calculated from the full scatter plots before binning. The solid black line is the 1 : 1 line and the 
dashed lines denote the expected uncertainty calculated from pre-launch analysis. 

Figure 10. Monthly mean aerosol optical thickness at 0.55 pm (heavy lines and left hand axes) and 
fraction of optical thickness contributed by fine mode at 0.55 pm ( t hin lines and right hand axes) for the 
year 2001 . Blue denotes MODIS ocean retrievals. Red denotes MODIS land retrievals, and black denotes 
AERONET. Four land stations are shown in (a) and four ocean stations are shown in (b). The MODIS 
values are calculated from Level 3 daily statistics and represent a 3 degree latitude by 3 degree longitude 
box centered on the AERONET station. 

Figure 11. MODIS retrieved aerosol particle effective radius over ocean plotted against AERONET 
retrievals of the same parameter. Only points with AERONET x 0 ' 44 > 0. 1 5 are plotted. AERONET values 

are daily averages for the date of the MODIS overpass. The blue line represents the linear regression 
through the points. The solid black line is the 1:1 line and the dashed lines represent ±0. 1 0 pm. 27 1 co- 
located points are shown. 62% of these points fall within the dashed lines. 

Figure 12. MODIS aerosol optical thickness at 0.55 pm, representing global aerosol distribution on August 
12, 2001. The two dimensional color bar describes both magnitude of optical thickness (along bottom axis) 
and fraction of optical thickness contributed by smaller fine mode particles (along vertical axis). Blue 
indicates low aerosol loading. Red indicates heavy loading of small particles such as pollution and smoke. 
The greener tones indicate a greater percentage of large particles such as desert dust and sea salt. The 
image was created from the MODIS daily 10-km resolution data after smoothing the raw data with 
Gaussian filters applied both in the temporal and spatial domains. 



Table 1. Size Distribution Parameters, and Single Scattering Albedo Used in the MODIS Look Up 
Table for the Land Algorithm. 
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Table 2: Refractive Indices, Median, Standard Deviation and Effective Radius for 
the aerosol models used in the MODIS Lookup Table for the ocean algorithm. 
Models 1-4 are fine modes and Models 5-9 are coarse modes. 
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Table 3 Values of Assymetry parameter and Single Scattering Albedo for the 9 ocean models of 
Table 2. 
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0.5755 

0.5117 

0.4478 

0.3221 

0.1773 

0.1048 

0.0622 

2 

0.6832 

0.6606 

0.6357 

0.5756 

0.4677 

0.3685 

0.2635 

3 

0.7354 

0.7183 

0.6991 

0.6510 

0.5590 

0.4715 

0.3711 

4 

0.7513 

0.7398 

0.7260 

0.6903 

0.6179 

0.5451 

0.4566 

5 

0.7450 

0.7369 

0.7328 

0.7316 

0.7330 

0.7411 

0.7282 

6 

0.7770 

0.7651 

0.7503 

0.7358 

0.7314 

0.7461 

0.7446 

7 

0.8035 

0.7912 

0.7738 

0.7506 

0.7335 

0.7443 

0.7461 

8 

0.7534 

0.7200 

0.6979 

0.6795 

0.7129 

0.7173 

0.7190 

9 

0.7801 

0.7462 

0.7234 

0.7065 

0.7220 

0.7176 

0.7151 

Single seal 

ttering albedo 

1 

0.9735 

0.9683 

0.9616 

0.9406 

0.8786 

0.5390 

0.4968 

2 

0.9782 

0.9772 

0.9757 

0.9704 

0.9554 

0.8158 

0.8209 

3 

0.9865 

0.9864 

0.9859 

0.9838 

0.9775 

0.9211 

0.9156 

4 

0.9861 

0.9865 

0.9865 

0.9855 

0.9819 

0.9401 

0.9404 

5 

0.9239 

0.9358 

0.9451 

0.9589 

0.9707 

0.9753 

0.9774 

6 

0.8911 

0.9026 

0.9178 

0.9377 

0.9576 

0.9676 

0.9733 

7 

0.8640 

0.8770 

0.8942 

0.9175 

0.9430 

0.9577 

0.9669 

8 

0.9013 

0.9674 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

9 

0.8669 

0.9530 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 



Table 4. Contents of MODIS Aerosol Level 2 hdf file (MOD04): Time and geometric 
information. 


Name of product 

Dimensic 

status 

Longitude 

2D 

diagnostic 

Latitude 

2D 


Scan_Start_Time 

2D 


Solar_Zenith 

2D 

diagnostic 

Solar_Azimuth 

2D 


Sensor_Zenith 

2D 


Sensor_Azimuth 

2D 


Scattering_Angle 

2D 


Cloud_Mask_QA 

2D 



Dimension: 2 dimensional arrays of 204 x 135 are indicated. If the array is 3 dimensional, the 
values of the 3rd dimension are given. The 3rd dimension 

Table 5. Contents of MODIS Aerosol Level 2 hdf file (MOD04): Global land and ocean products, 
at 550 nm. 


Name of product 

Dimension 

status 

Optical_Depth_Land_And 

__Ocean 

2D 

validated 

OpticalJDepth_Ratio__Small 
_ Land_And_Ocean 

2D 

Not yet 
validated 

■BUM 

2D 

derived 
























Table 6. Contents of MODIS Aerosol Level 2 hdf file (MOD04): Land products. 


Name of product 

Dimension 

Status 

Corrected_Optical_Depth 

JLand 

0.47, 0.55, 
0.66 fiin 

Validated 

Optical_Depth_Ratio 

_Small_Land 

0.55 pm 

Not yet valid 

Mass_Concentration„Land 

2D 

Derived 

Angstrom_Exponent_Land 

0.66/0.47 

Not yet valid 

Reflected_FluxJLand 

0.47, 0.55, 
0.66 um. 6 

Derived 

T ransmittedjdux JLand 

0.47, 0.66 

|rm 

Derived 

AerosolJTypeJLand 

2D 

Diagnostic 

Continental_Optical_Depth 

JLand 

0.47, 0.66 
pm 

Diagnostic 

EstimatedJJncertainty_Land 

0.47, 0.66 
pm 

Diagnostic 

Mean_Reflectance_Land_All 

0.47, 0.66, 
2.13 pm 

Diagnostic 

S tandard JDeviati on 

__Reflectance_Laiid_All 

0.47, 0.66, 
2.13 pm 

Diagnostic 

CloudJFractionJLand 

2D 

Diagnostic 

Number_Pixels_Percentile 

_Land 

0.47, 0.66 
pm 

Diagnostic 

Mean_Reflectance_Land 

0.47, 0.66, 
0.87,2.13, 
3.75 pm 

Diagnostic 

STD_Reflectance__Land 

0.47, 0.66, 
0.87,2.13, 
3.75pm 

Diagnostic 

Quality_Assurance_Land 

See QA 
plan 

Diagnostic 

Path_Radiance_Land 

0.47, 0.66 
pm 

Experimenta] 

Error JP ath_Radi anc e„Land 

0.47, 0.66 
pm 

Diagnostic 

Critical_Reflectance_Land 

0.47,0.66 

pm 

Experimenta! 

Error_CriticaLReflectance 

_Land 

0.47,0.66 

pm 

Diagnostic 

QualityWeightJPath_Radianc( 

0.47,0.66 

pm 

Experimenta! 

QualityWeight_Critical 

_Reflectance_Land 

0.47, 0.66 
pm 

Experimenta] 

Quality_Assurance_Crit_Ref 

_Land 

0.47, 0.66, 
0.87,2.13, 
3.75 pm 

Diagnostic 









Table 7. Contents of MODIS Aerosol Level 2 hdf file (MOD04): Ocean products. 


Name of product 

Dimension 

Status 

Effective_Optical_Depth_Average_Ocean 


Validated 

OpticalJDepth__SmalLAverage 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63, 2.13 jim 

Not yet validated 

Optical_D epth__Large_A verage 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63,2.13 jim 

Not yet validated 

Effective_Radius_Ocean 

Best, Average 

Validated* 

Optical_D epth_Ratio_Small_Ocean_0 . 8 6micron 

Best, Average 


Mass_Concentration_Ocean 

Best, Average 

Derived 

Cloud_Condensation_Nuclei_Ocean 

Best, Average 

Derived 

Angstrom_Exponent_l_Ocean 

0.55/0.87 


Angstrom_Exponent_2_Ocean 

0.87/2.13 

Not yet validated 

Reflected_Flux_Average_Ocean 


Derived 

Transmitted__Flux_Average_Ocean 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63,2.13 jun 

Derived 

A symmetry JF actor_Average_Ocean 

0.47, 0.55, 0.66, 0.87 * 
1.24, 1.63, 2.13 pm 

Derived 

B ackscattering_Ratio_Average_Ocean 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63, 2.13 jim 

Derived 

Solution Jndex_Ocean_Small 

Best, Average 


Solution„Index_Ocean_Large 

Best, Average 


Least_Squares_Error_Ocean 

Best, Average 


OpticalJDepth_by__models_ocean 

9 models 


Effective_0ptical_Depth_Best_0cean 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63, 2.13 Jim 

Diagnostic 

Optical__Depth_Sniall_B est 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63, 2.13 jim 

Diagnostic 

OpticalJDepthJLargeJBest 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63,2.13 Jim 

Diagnostic 

Reflected_Flux_Best_Ocean 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63,2.13 Jim 

Diagnostic 

Transmitted_Flux_Best_Ocean 

0.47, 0.55, 0.66, 6.87 
1.24, 1.63,2.13 Jim 

Diagnostic 

Asymmetry_Factor_Best_Ocean 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63, 2.13 jim 

Diagnostic 

Backscattering_Ratio_Best_Ocean 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63,2.13 Jim 

Diagnostic 

Cloud_Fraction_Ocean 

2 dimensional 

Diagnostic 

Number_Pixels_Used_Ocean 

2 dimensional 


Mean_Reflectance„Ocean 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63, 2.13 jim 

Diagnostic 

STD_Reflectance_Ocean 

0.47, 0.55, 0.66, 0.87 
1.24, 1.63, 2.13 jim 

Diagnostic 

QuaHty_Assurance__Ocean 

See QA plan 

Diagnostic 




















































Table 8. Number of retrievals (N), percentage of retrievals (%) falling within expected uncertainty 
for each derived wavelength, average optical thickness from AERONET (x°' 55 ) and percent 
difference between MODIS and AERONET at 0.55 pm (Diff) for the global data set and for each 
region, land and ocean separately. 


Region 

N 

% 0.47 

% 0.55 

% 0.66 

% 0.87 

x 0 - 55 

Diff 

All land 

5906 

61 

68 

71 


0.18 

41 

China_land 

205 

71 

75 

76 


0.28 

17 

India_land 

70 

70 

90 

79 


0.46 

-2 

Mediterranean^ and 

781 

65 

67 

66 


0.20 

32 

Alaska/ Canada_land 

178 

49 

57 

67 


0.10 

122 

South Afiica_land 

r"347 

79 

86 

84 


0.19 

-10 

South America_land 

762 

63 

72 

75 


0.18 

21 

West USJand 

901 

56 

61 

63 


0.14 

53 

East US _land 

1385 

55 

64 

68 


0.17 

54 

East Europe_land 

192 

79 

83 

72 


0.24 

-10 

West Europe_land 

793 

66 

72 

77 


0.17 

48 


Region 


N 


% 0.47 % 0.55 % 0.66 % 0.87 


0.55 


Diff (%) 


All ocean 

Indian_ocean 

Asian Pacific_ocean 
Pacific Island_ocean 
West Medit_ocean 
East Medit_ocean 
Sahara n_ocean 
Atlantic Isles_ocean 
Australia_ocean 
North Europe_ocean 
Caribbean_ocean 
East Pacific_ocean 
US Atlantic_ocean 


2052 

47 

57 

163 

334 

205 

184 

146 

70 

150 


62 

64 

56 
70 
52 

57 
58 . 

64 
83_ 

65 
62 
52 
72 


66 

77 

53 

74 

62 

63 

56 

7J_ 

8j_ 

72 

67 
61_ 

68 


70 

77 

60 

79 

68 

II 

51_ 

II 

83 

8J_ 

68 

69 

70 


0.18 

0.16 

0.21 

0.08 

0.21 

0.23 

0.31 

0.13 

0.05 

0.16 

0.14 

0.18 

0.15 


1 _ 

13 

- 6 _ 

- 6 _ 

£7_ 

_1_ 

8 _ 

2 _ 

A 

20 

- 6 _ 

1 


242 

160 

288 

























MODIS Over Land Algorithm 

All procedures applied to individual boxes of 20 x 20 pixels at 500 m resolut ion ( 1 0 km at nadir) ^ 

Identify and mask all water, cloudy and snow/ice pixels. 

Separate all pixels that are outside of 0.0 1< p 213 ^ 0.25 
Discard brightest 50% and darkest 20% of pixels defined with p 0 66 




Output x, Flux and other parameters 


* Pure dust 

Use dust LUT to make retrieval 
Only IF 0.15 <p 2 . 13 < 0.25 


Out put T, Flux and other parameters 


Make dust and non-dust 
Retrievals and combine 
By weighting with 

/-» 0.66 / 0.47 

Ko ' ro 


Mixed 


0.72 < AND <0.90 f( 


Path A - Set QC=3 

For all identified ‘dark’ pixels, 
Calculate mean p 0 47 , p 066 
p s 0 - 47 =0.25 p 2 I3 ; p s 0 - 66 =0.50 p 213 


Continental model LUT 


t 0 .47 and x 066 


Calculate path radiance (p 0 ) at 0.47 and 0.66 
Spectral dependence of path radiance 
distinguishes dust from non-dust. 


Path B - Set QC=0 

Now collect pixels with p 2 13 > 0.25 
previously set aside. Count number ot 
pixels (N) with 

0.0 1< p 2 13 < 0.25 f(p,|J 0 ) 

AND p 2 13 < 0.40 

^ "X. 

yes 


fill values 


p s oj7 =0.25 p 213 ; p s 0 66 =0.50 p 


2.13 


Continental model LUT 


t 0.47 


Figure Flowchart illustrating the derivation of aerosol over land. 











January 2003 ■■■^■e==shb April 2003 



July 2002 October 2002 


Figure 2. Monthly mean plots of fraction of total aerosol optical thickness attributed to 'non 
dust' or fine mode aerosol over land. Fraction 1.0 indicates all fine mode. Fraction of 0.0 
indicates all coarse mode. 


Urban/Industrial(0.96) 




V 




B," / H 


smoKe - moderate ^ 
absorption (0.90) — 


Hig 

Sea 


Figure 3. Distribution of the non-dust models 
Single scattering albedo values given in parent 
aerosol model. Solid black regions use the me 
checkerboard pattern in east Asia denotes a rej 
Africa is divided by region and season. North 
(November-May) the strong absorption aeroso 
model is used the remainder of the year. South 
October when the strong absorption model is u 
algorithm uses the urban/industrial model. 





Year 2000 Day 233 Time 0835 
Location: southern Africa, east coast 



Figure 4. MODIS-derived aerosol optical thickness at 0.55 pm for an image of the east coast of 
southern Africa. Top panel uses the traditional dark target method described by Path A in 
Figure 2. The bottom panel shows the results after extending the retrieval to brighter surfaces 
as described by Path B of Figure 2. By extending to brighter surfaces the number of retrievals 
over land in this image increases from 7060 to 17,849. 






MODIS Aerosol Over Ocean Algorithm 

All procedures applied to individual boxes of 20 x 20 pixels at 500 m resolution (10 km at nadir) 


Entire box clear of land pixels according to MOD35 lkm mask' No? 


^es? ► Continue masking pixel by pixel 

( 1 ) spatial variability : stdev of p ®' 55 calculate 
from 3x3 centered on pixel >0.0025 then 
all 9 pixels labeled ‘cloudy’. 

(2) dust call back: if p 0 47 / p 0 66 < 0.75, then 
dust, use even if variability is high. 

(3) if p 0,47 > 0.40 then ‘cloudy’ 

(4) IR tests : if any of 3 specific MOD35 tests 
indicate ‘cloudy’, then ‘cloudy’ 

(5) 1.38 cirrus tests : cloudy if p 1 38 / p * 24 > q 3 
OR if 0.01 < p 1 38 / p 1 24 ^ 0.30 

AND p 138 > 0.03 
but, if 0. 10 £ p 1 - 38 / p 1 24 £ 0.30 

AND 0.01 < p 1 38 < 0.03 then 
‘not cloudy’ but QC=0 

( 6 ) sediment mask: 



Discard brightest 25% and 
darkest 25% defined with p ^ 8 ^ 


'Flu values^ Are there 

at least 10 
./good’ pixels^ 


EXIT 


yes, 


Calculate mean and stdev of p in 7 A,’s 




Heavy dust 
over glint. 
QC=0 



consistent 


EXIT 


i 

Do Inversion 

For all 20 combinations 

of 1 fine mode and 1 coarse mode. 

Find T) and x such that 

p m ^ - p tot ^ is minimum over all X’s 
where p m is MODIS measured reflectance 
and p to j = r\ + ( l-r|) p c ^ from the LUT. 
Choose which of the 20 combinations 
has least error (e) 

Note in minimizing e, wavelengths are 
weighted for information content and QC. 

Chosen LUTs ► x, Flux in 7 A.’s, r_eff etc. 

Best solution is least error. 

Average solution is either: 

(a) average of all solutions where e < 3% 

(b) if no e < 3%, then average of best 3 

Final Checking 


If T ^ 5.0 (putput x, flux in 7 X’s, r_eff etc. 


Except For Dust 
when x < 0.7 (Till value 


Ifx>5.0 For Dust, setx=5.0 




Output x=5.0. Flux in 7 A,’s, 
r_eff QC = 0 


For Not Dust 




F-ijore. 5~ 









Figure 6. Examples of MODIS aerosol products at the Level 2 stage (MOD04). The data 
represent a 5-minute granule collected on July 7, 2002 from 1835-1840 UTC. The upper left 
panel is a true color image created from Level lb reflectances. The upper right panel is aerosol 
optical thickness at 550 nm, The lower left panel is reflected flux at 550 nm. The lower right 
panel is fraction of fine mode. 





Figure 7. The distribution of the 132 AERONET stations used to validate MODIS land and 
ocean aerosol retrieval algorithms. 





y = 0.068 + 0.78 x R= 0.80 



y = 0.059 + 0.70 x R = 0.68 



Figure 8. MODIS aerosol optical thickness retrievals over land at 470 nm (blue), 550 nm (green 
) and 660 nm (red) as a function of AERONET observations co-located in space and time. The 
data were sorted according to AERONET aerosol optical thickness and averaged for every 300 
points. At higher optical thickness where the data beome sparser, fewer points are used in the 
average, as indicated. The standard deviation in each bin is shown by error bars. The 
regression equations given at the top of each plot were calculated from the full scatter plots 
before binning. The solid black line is th 1:1 line and the dashed lines denote the expected 
uncertainty calculated from pre-launch analysis. 




R = 0.92 


y = 0.008 + 0.95 x 



y = 0.005 + 0.94 x R = 0.91 



y = 0.016 + 0.87 x R= 0.91 



Figure 9. MODIS aerosol optical thickness retrievals over ocean at 550 nm (green ) and 660 nm 
(red) and 870 nm (black) as a function of AERONET observations co-located in space and time. 
The data were sorted according to AERONET aerosol optical thickness and averaged for every 
100 points. At higher optical thickness where the data become sparser, fewer points are used in 
the average, as indicated. The standard deviation in each bin is shown by error bars. The 
regression equations given at the top of each plot were calculated from the full scatter plots 
before binning. The solid black line is the 1:1 line and the dashed lines denote the expected 
uncertainty calculated from pre-launch analysis. 










Figure 10. Monthly mean aerosol optical thickness at 0.55 pm (heavy lines and left hand axes) 
and fraction of optical thickness contributed by fine mode at 0.55 pm ( thin lines and right hand 
axes) for the year 2001. Blue denotes MODIS ocean retrievals. Red denotes MODIS land 
retrievals, and black denotes AERONET. Four land stations are shown. The MODIS values are 
calculated from Level 3 daily statistics and represent a 3 degree latitude by 3 degree longitude 
box centered on the AERONET station. 


fine mode AOT/total AOT fine mode AOT/tota! AOT 



Aerosol Optical Thickness Aerosol Optical Thickness 




Figure 10b. 


fine mode AOT/total AOT fine mode AOT/total AOT 
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Effective Radius Comparison 
t 044 >0.15 



AERONET Effective Radius ( pm) 


Figure 11. MODIS retrieved aerosol particle effective radius over ocean plotted against 
AERONET retrievals of the same parameter. Only points with AERONET t 044 > 0.15 are 
plotted. AERONET values are daily averages for the date of the MODIS overpass. The blue 
line represents the linear regression through the points. The solid black line is the 1:1 line and 
the dashed lines represent ±0.10 pm. 271 co-located points are shown. 62% of these points fall 
within the dashed lines. 





Figure 12. MODIS aerosol optical thickness at 550 nm, representing global aerosol distribution 
on August 12, 2001. The two dimensional color bar describes both magnitude of optical 
thickness (along bottom axis) and fraction of optical thickness contributed by smaller fine mode 
particles (along vertical axis). Blue indicates low aerosol loading. Red indicates heavy loading 
of small particles such as pollution and smoke. The greener tones indicate a greater percentage 
of large particles such as desert dust and sea salt. The image was created from the MODIS daily 
10-km resolution data after smoothing the raw data with Gaussian filters applied both in the 
temporal and spatial domains. 





